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Abstract 

Fuel  sensitivity  experiments  form  an  important  feature  of  testing  solid  oxide  fuel  cell  systems.  In  this  paper  we  describe  a  fuel  sensitivity 
test  and  its  relation  to  fuel  utilization.  The  sensitivity  of  the  measured  terminal  voltage  of  a  tubular  solid  oxide  fuel  cell  to  the  fuel  utilization 
provides  important  information  about  ‘leaks’  through  the  cell.  Such  leaks  could  arise  from  pinholes  in  the  electrolyte  or  from  other  sources. 
First  a  simple  analytical  model  is  presented  which  is  then  refined  using  a  numerical  simulation.  Such  tests  also  provide  a  methodology  to 
quantify  the  leaks  present  in  the  cell  which  can  ultimately  have  an  effect  on  fuel  cell  efficiency. 
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1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  are  being  widely  in¬ 
vestigated  for  application  in  environmentally  clean  power 
generation.  SOFCs  have  many  important  advantages  over 
conventional  power  generation  systems  including  high  ef¬ 
ficiency,  noiseless  operation,  modular  construction  and 
essentially  zero  NO*  and  SO*  emissions  [1,2].  Due  to  their 
high  efficiencies  they  also  emit  a  much  smaller  amount  of 
the  greenhouse  gas  CO2  compared  to  conventional  power 
generation  systems.  Seal-less  tubular  SOFCs  are  at  the 
forefront  of  this  technology.  In  this  paper  we  present  an 
experimental  and  analytical  study  of  the  sensitivity  of  the 
tubular  SOFC  voltage  to  the  fuel  utilization.  The  cell  voltage 
at  constant  current  density  has  been  measured  as  a  func¬ 
tion  of  fuel  utilization.  It  is  demonstrated  that  the  observed 
behavior  can  be  explained  using  a  simple  analytical  model. 


2.  Fuel  utilization  sensitivity  experiments 

Fuel  utilization  sensitivity  tests  are  routinely  performed 
on  tubular  SOFCs.  Such  tests  are  used  to  make  a  qualitative 
diagnostic  determination  of  the  operational  fuel  utilization 
of  single  cells.  In  a  fuel  utilization  sensitivity  experiment, 
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the  cell  current  (or  current  density)  is  held  constant  and  the 
change  in  cell  voltage  is  measured  as  a  function  of  the  fuel 
utilization  over  a  wide  range.  When  the  current  density  is 
held  constant,  the  total  fuel  flow  to  the  cell  is  not  constant 
but  is  inversely  proportional  to  the  fuel  utilization.  The  fuel 
utilization  is  defined  as  the  fraction  of  the  input  fuel  which  is 
electrochemically  utilized,  i.e.  the  fuel  utilization  U  is  given 
by  the  following  equation: 


In  the  above  equation,  Qqc  (mol/s)  is  the  part  of  the  fuel  that 
is  electrochemically  utilized  in  generating  electricity  and  Qt 
(mol/s)  the  total  fuel  input  to  the  cell.  The  electrochemically 
utilized  fuel  can  be  related  to  the  cell  current  using  Faraday’s 
law  as  follows: 


I  _J  A 
IF  ~  2F 


where  I  is  the  cell  current  (A),  J  the  average  current  density 
(A/cm2),  A  the  electrode  area  (cm2),  and  F  the  Faraday’s 
constant  (C/mol).  It  is  important  to  note  that  in  general  the 
current  density  /  in  a  tubular  fuel  cell  varies  as  a  function  of 
axial  position  and  the  current  density  in  Eq.  (2)  J  represents  a 
position  averaged  current  density.  Fig.  1  shows  the  geometry 
of  a  single  cell  test  article.  The  terminal  voltage  of  a  tubular 
SOFC  is  given  by  the  following  equation: 

Fterm  —  FNernst(x)  RJ(x)  7j(v)  (3) 
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Nomenclature 

A 

fuel  electrode  (anode)  area  (cm2) 

AE/EL 

at  the  cathode  (air  electrode)-electrolyte 
interface 

Dlu 

effective  H2-H2O  diffusivity  in  the 
anode  (cm2/s) 

°eff 

effective  N2-O2  diffusivity  in  the 
cathode  (cm2/s) 

EL/FE 

at  the  electrolyte-anode  (fuel  electrode) 
interface 

F 

Faraday  constant  (96  487  C) 

I 

cell  current  (A) 

J 

current  density  (mA/cm2) 

Ki 

equilibrium  constant  for  H2-H2O  reaction 
(atm-0-5) 

l 

total  electrochemically  active  length  of 
tubular  SOFC  (cm) 

air 

P02 

cathode  (air  electrode)  side  oxygen  partial 
pressure  (atm) 

„fuel 

Po2 

anode  (fuel  electrode)  side  oxygen  partial 
pressure  (atm) 

Q 

total  fuel  flow  (mol/s) 

Q&c 

fuel  flow  which  is  electrochemically 
utilized  (mol/s) 

r 

radial  position  within  the  cell  measured 
from  axis  of  the  tubular  SOFC  (cm) 

R 

universal  gas  constant  (8.314  J/(molK)) 

Ri 

area  specific  resistance  of  the  tubular 

SOFC  (fi  cm2) 

G 

anode  thickness  (cm) 

tc 

cathode  thickness  (cm) 

u 

fuel  utilization  (0-1) 

^Nernst 

local  Nernst  potential  (mV) 

Verm 

cell  terminal  voltage  (mV) 

axial  position  measured  from  the  fuel  inlet 
point  (closed-end)  of  the  tubular  SOFC  (cm) 

Greek  letter 

T] 

electrode  concentration  polarization  (mV) 

In  the  above  equation,  Vterm  is  the  terminal  voltage  which  is 
constant  and  independent  of  axial  position  at  a  given  current 
density  and  fuel  utilization,  VNernst  the  local  Nernst  poten¬ 
tial  at  axial  position  x,  R  the  area  specific  resistance  (£2  cm2) 
which  is  independent  of  position,  J  the  current  density  which 
is  a  function  of  cell  axial  position  and  r)  (v)  the  sum  of  var¬ 
ious  electrode  polarization  losses  at  the  cathode  and  anode 
which  is  also  expected  to  be  a  function  of  cell  axial  position. 
It  is  important  to  note  that  the  terminal  voltage  is  indepen¬ 
dent  of  cell  axial  position  due  to  the  presence  of  highly  con¬ 
ductive  axially  positioned  current  collecting  bus-bars  which 
act  as  equipotential  surfaces.  Taking  averages  of  every  term 
in  Eq.  (3)  the  terminal  voltage  can  be  written  in  terms  of  the 


Post-combustion 


average  values  as  follows: 

Verm  —  ^Nernst  JR  0  (4) 


The  average  Nernst  potential  is  the  position  averaged  Nernst 
potential.  In  other  words 

1  fx=l 

^Nernst  —  T  /  VNernstCO  d-T  (5) 

1  Jx= 0 


Such  an  average  is  necessary  since  the  Nernst  potential  at 
a  certain  axial  location  along  the  cell  length,  under  the  as¬ 
sumption  that  the  oxygen  utilization  is  negligible,  i.e.  air  is 
fed  to  the  cell  far  in  excess  of  the  stoichiometrically  required 
amount  is  given  by 


^Nernst  00  — 


4  F 


In  the  above  equation,  the  partial  pressure  of  oxygen  on  the 
airside  p ^  is  0.21  atm  and  is  independent  of  axial  position 
due  to  negligible  oxygen  utilization.  On  the  fuel  side  the  oxy¬ 
gen  partial  pressure  p^(x)  is  determined  by  the  H2-H2O 
equilibrium  and  is  a  function  of  cell  axial  position.  The  oxy¬ 
gen  partial  pressure  in  the  fuel  can  then  be  related  to 
the  partial  pressure  of  water  vapor  and  hydrogen  on  the  fuel 
side  p^o  M  and  Puf  M*  from  the  H2-H2O  equilibrium 
through  the  following  relationship: 


where  Afz  is  the  equilibrium  constant  for  the  H2-H2O.  Hence¬ 
forth,  the  superscript  fuel  will  be  dropped  from  p^0  (v)  and 
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p^e\x)  as  these  species  are  present  only  on  the  fuel  side  of 
theSOFC. 

Using  Eq.  (7)  one  can  write  for  the  oxygen  partial  pres¬ 
sures  in  the  fuel  stream  at  the  fuel  inlet  and  exit  points  as 


itly  addressed  in  the  numerical  simulations  presented  later 
in  the  paper.  Under  the  foregoing  assumptions 


dVterm  _  d  V^sfernst 

d  U  ~  d  U 


(14) 


=  0)]  =  2  In 


Ph2o(x  =  Q) 

_  pu2(x  =  0)  _ 


2  In  (Kt) 


PH  2o(x  =  L) 
Ph2(x  =  L) 


2  In  (Ki) 


(8) 

(9) 


We  now  make  an  assumption  that  the  logarithm  of  oxygen 
partial  pressure  on  the  fuel  side  is  a  linear  function  of  ax¬ 
ial  position.  This  assumption  which  is  accurate  for  high  av¬ 
erage  current  densities  is  based  on  a  numerical  simulation 
described  later  in  the  paper  which  predicts  that  the  Nernst 
potential  is  linear  with  cell  axial  position  (see  Fig.  3).  Under 
this  assumption,  the  mean  In (pfQ^1)  across  the  length  of  the 

cell  is  simply  given  by  the  arithmetic  average  of  \n[pf^\x  = 

0)]  and  ln[/?Q^  (x  =  L)].  Using  Eqs.  (8)  and  (9)  the  average 
oxygen  partial  pressure  at  high  average  current  densities  on 
the  fuel  side  is  given  by 


Ph2o(x  =  0) 
_  pn2  (x  =  0) 


+  In 


Ph2o(x  =  L) 
pu2(x  =  L) 


-2  In  (Ki) 


(10) 


The  average  Nernst  potential  at  high  current  densities  is  then 
given  by 


VNernst,h 


-4#') 

4f  \p“/ 


(ID 


where  Pq^\  is  given  in  Eq.  (10),  where  the  subscript  h  stands 
for  high  average  current  density.  In  the  above  equation  we 
have  assumed  that  the  air  side  oxygen  partial  pressure  does 
not  change  with  axial  position  since  the  air  supplied  p™n  to 
the  cell  during  these  tests  is  several  times  (10-15  times)  the 
amount  needed  for  complete  stoichiometric  oxidation  of  the 
fuel. 

The  water  vapor  and  hydrogen  partial  pressures  in  the  exit 
of  the  fuel  cell  can  be  given  by 


Ph2o(x  =  L)  =  Ph2o(x  =  0)  +  Upu2(x  =  0)  (12) 


pu2(x  =  L)  =  pu2(x  =  0) [1  -  U] 


(13) 


where  U  is  the  fuel  utilization. 

We  now  assume  that  the  second  and  third  terms  in  Eq.  (4) 
are  independent  of  fuel  utilization.  Since  the  tests  are  per¬ 
formed  at  constant  average  current  density  and  since  the 
ohmic  resistance  of  the  cell  is  not  expected  to  change  with 
fuel  utilization  neglecting  the  second  term  is  well  justi¬ 
fied.  However,  the  third  term  may  well  depend  on  the  fuel 
utilization.  Thus  this  analysis  can  be  regarded  as  a  zeroth 
order  approximation.  The  effect  of  fuel  utilization  on  the 
third  term,  which  represents  polarization  losses,  is  explic- 


Substituting  Eqs.  (10),  (12)  and  (13)  in  Eq.  (11)  and  dif¬ 
ferentiating  the  Nernst  potential  with  respect  to  the  fuel 
utilization  U,  and  using  Eq.  (14)  we  get 

dVterm, h  _  dVNernst5h 

d  U  ~  d U 

_  RT  (  1  \ 

4 F  \(l-U)[pu2o(x=0)  +  Upu2(x  =  0)] / 

(15) 


At  low  average  current  densities,  a  significant  length  of  the 
cell  would  be  operating  close  to  the  exit  Nernst  potential.  In 
other  words,  most  of  the  fuel  will  be  consumed  very  close 
to  the  fuel  inlet,  which  is  the  closed-end  of  the  cell  and  the 
Nernst  potential  would  decrease  sharply  from  the  inlet  value 
to  the  exit  value  within  a  short  length  from  the  inlet.  This 
can  indeed  be  seen  in  Fig.  3  for  an  average  current  density 
value  of  lOmA/cm2.  Thus,  the  exit  Nernst  potential,  which 
is  strictly  a  function  of  fuel  utilization  at  constant  temper¬ 
ature,  would  be  attained  very  close  to  the  closed-end  of  the 
cell.  Because  a  very  large  portion  of  the  cell  experiences 
a  Nernst  potential  close  to  that  of  the  exit,  at  a  very  low 
average  current  density,  the  average  Nernst  potential  can 
be  approximated  to  that  of  the  exit  Nernst  potential.  In  this 
case  the  average  oxygen  partial  pressure  across  the  cell, 
Po2\  (where  the  subscript  1  denotes  low  current  density)  is 
given  by 

ln(/4el,)  =  2  In 


PH2o(x  =  L ) 

_  Ph2  ( x  =  L) 


-2  In  (Ki) 


(16) 


The  average  Nernst  potential  in  this  case  is  given  by 


VNernst,! 


4  F 


(17) 


Substituting  Eqs.  (12),  (13)  and  (16)  in  Eq.  (17)  and  using 
Eq.  (14)  we  get 


dVterm,! 
d  U 


d  VjSfernst,  1 

d U 

RT  /  1  \ 

IF  V(1  -  U)[pn2o(x  =  0)+Upn2(x  =  0)] / 

(18) 


Comparison  of  Eqs.  (15)  and  (18)  reveals  that  the  derivative 
of  the  terminal  voltage  with  respect  to  fuel  utilization  at 
low  average  current  density  is  twice  that  at  high  average 
current  densities. 

Fig.  2  shows  the  plots  of  the  absolute  values  of  dVterm/d  U 
versus  U  over  a  range  of  fuel  utilization  values  for  both  low 
average  current  density  and  high  average  current  density 
cases  along  with  experimentally  measured  data  for  average 
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Fig.  2.  Fuel  utilization  sensitivity  curves  for  average  current  densities  of  100  and  500  mA/cm2  and  analytical  model  predictions. 


current  densities  of  100  and  5 00 mA/cm2.  As  can  be  seen 
from  Fig.  2,  although  the  experimental  data  do  not  match 
the  model  developed  here  perfectly,  the  overall  trend  has 
good  qualitative  agreement  with  observed  data.  Clearly  the 
model  needs  to  be  refined  to  accommodate  the  effect  of 
ohmic  losses  and  polarization  losses  to  get  a  more  accurate 
fit  to  the  experimental  observations.  This  is  addressed  in 
what  follows. 

Here  we  present  a  refinement  to  the  above  model  which 
includes  the  effect  of  the  second  and  third  terms  in  Eqs.  (3) 
and  (4).  Further  we  show  that  this  modification  also  allows 
one  to  use  the  fuel  sensitivity  experiment  as  a  tool  to  estimate 
leaks  in  tubular  SOFCs. 

The  polarization  loss  rj(x)  is  a  function  of  position  and  is 
the  sum  of  polarization  losses  on  the  air  electrode  and  fuel 
electrode  sides,  i.e. 

r)(x)  =  rfr{x)  +  )7fuel(x)  (19) 


The  air  electrode  and  fuel  electrode  side  polarization  losses 
can  in  turn  be  written  as: 


•  RT 
rj  (jc)  =  —  In 
1  4  F 


air 
P  Oo 


AE/EL.  . 

IPo ,  (*) 


(20a) 


)?fuel(x) 


(20b) 


The  experimental  data  presented  here  have  been  mea¬ 
sured  at  1000  °C.  At  this  temperature,  we  have  found  the 
charge-transfer  polarization  losses  in  our  cells  are  negli¬ 
gible.  Thus  the  polarization  losses  in  Eq.  (19)  are  limited 


to  concentration  polarization  losses  as  given  in  Eqs.  (20a) 
and  (20b).  The  oxygen  partial  pressures  at  the  elec¬ 
trode/electrolyte  interfaces  p^^EL  (x)  and  /?q^FE  (jc)  can  be 
related  to  the  electrode  dimensions  and  effective  binary  gas 
diffusion  coefficients  as  outlined  by  Kim  et  al.  [3]  and  Chan 
et  al.  [4].  The  calculations  of  Kim  et  al.  [3]  and  Chan  et  al. 
[4]  pertain  to  a  planar  SOFC  with  negligible  fuel  utiliza¬ 
tion.  However,  the  calculations  in  the  present  work  are  for 
the  case  of  a  tubular  SOFC  with  a  large  electrode  area  with 
practical  fuel  utilization.  The  oxygen  partial  pressure  at  the 
cathode-electrolyte  interface  is  determined  by  relating  the 
oxygen  flux  Jo2  (mol/(cm2  s))  through  the  electrolyte  to  the 
current  density  J  (A/cm2),  i.e. 


^eff  ^7*02 
RT  dr 


(21) 


It  is  important  to  note  that  in  Eq.  (21),  the  current  density 
(and  thus  the  oxygen  flux  through  the  electrolyte),  and  the 
radial  oxygen  partial  pressure  gradient  are  a  function  of  cell 
axial  position.  Linearizing  Eq.  (21)  we  can  write 

AE/EL 

>2  P02 


4  F 

J(x)  =  -  —  Dc 


air 

P0: 


(x) 


RT 


eff 


tc 


(22) 


In  Eq.  (22)  p ^  is  the  oxygen  partial  pressure  in  the  bulk 
air  stream,  Dceff  the  effective  binary  diffusivity  of  N2-O2  in 
the  gas  phase  and  tc  the  air  electrode  thickness.  Recognizing 
that  the  current  density  is  a  function  of  cell  axial  position, 
the  interfacial  oxygen  partial  pressure  can  be  obtained  from 
Eq.  (22)  as  follows: 


EL/AE 

Po2 


M  =  Pol 


RTtc 

4FDC 

^rLJcff 


I /Ml 


(23) 
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Similarly  on  the  fuel  side 

pf /EL  (x)  =  PU2(x)  -  oF/ja  I I  (24a) 

^eff 

PFfoL  (*)  =  PHioix)  +  1  m  |  (24b) 

These  can  then  be  substituted  into  Eqs.  (20a)  and  (20b)  and 
the  results  of  the  substitution  into  Eq.  (3).  This  gives  the 
terminal  voltage  of  the  cell  as  a  function  of  the  local  cur¬ 
rent  density.  As  defined  earlier,  the  average  current  density 
J  is  defined  as  the  total  cell  current  divided  by  the  anode 
area,  i.e.  1/ A.  For  a  given  temperature,  fuel  utilization  and 
average  current  density,  one  can  iterate  to  obtain  the  current 
density  distribution  and  the  cell  terminal  voltage.  In  this  cal¬ 
culation  the  cell  is  divided  into  a  given  number  of  elements 
(which  can  be  varied),  and  at  a  given  average  current  den¬ 
sity,  for  each  of  these  elements,  a  local  fuel  utilization  and 
an  elemental  current  can  be  defined.  The  local  current  den¬ 
sity  and  the  local  interfacial  partial  pressures  of  the  various 
gases  can  be  calculated  as  a  function  of  cell  axial  position 
v  through  Eqs.  (1 1)— (13).  These  results  can  then  be  substi¬ 
tuted  in  Eqs.  (6)-(8)  to  calculate  the  Nernst  potential  and 
the  concentration  polarization  voltage  losses  on  the  cathode 
and  anode  sides  as  a  function  of  cell  axial  position  v.  In  the 
calculation,  R  which  is  the  area  specific  resistance  (£2  cm2) 
of  the  cell  and  the  diffusivities  £Eff  and  D^ff  can  be  var¬ 
ied  to  obtain  a  good  fit  with  the  experimentally  observed 
voltage-current  density  curves  (not  shown  here).  Further, 
the  fuel  leak  from  a  tubular  SOFC  can  be  incorporated  in 


the  analysis  by  including  a  subtractive  term  as  a  variable  pa¬ 
rameter  in  the  current  density  in  Eqs.  (23),  (24a)  and  (24b). 
This  subtractive  term  is  directly  proportional  to  the  magni¬ 
tude  of  fuel  leak. 

The  Nernst  potential  as  a  function  of  cell  axial  position 
calculated  using  the  numerical  scheme  described  above  is 
shown  in  Fig.  3  for  two  current  densities,  lOmA/cm2  (low) 
and  520  mA/cm2  (high).  The  dotted  line  for  the  520  mA/cm2 
case  represents  the  linear  approximation  of  the  variation  of 
Nernst  potential  with  respect  to  cell  axial  position  for  the 
high  average  current  density  case. 

As  can  be  seen  from  Fig.  3,  the  assumption  of  a  linear 
variation  in  the  Nernst  potential  as  a  function  of  axial  posi¬ 
tion  in  the  high  average  current  density  case  is  well  justified. 
Fig.  3  shows  the  results  of  the  variation  of  Nernst  poten¬ 
tial  at  520  mA/cm2  obtained  from  the  numerical  simulation 
and  compares  it  to  a  straight  line  fit.  Further,  the  assump¬ 
tion  of  an  average  Nernst  potential  (and  therefore  logarithm 
of  the  fuel  side  oxygen  partial  pressure)  equal  to  the  exit 
Nernst  potential  for  the  low  average  current  density  case  is 
also  equally  well  justified.  A  subtle  point  needs  to  be  made 
here  regarding  the  values  of  the  inlet  Nernst  potential  in  the 
case  of  high  average  current  density  and  low  average  current 
density  cases.  From  Fig.  3,  it  is  seen  that  the  inlet  Nernst  po¬ 
tential  is  ^980  mV  for  the  high  average  current  density  case 
which  is  close  to  what  should  be  expected  for  fuel  of  com¬ 
position  89%  H2  and  11%  H2O.  However,  for  the  low  aver¬ 
age  current  density  case  it  is  found  to  be  ^920  mV.  At  first 
sight,  this  is  puzzling  as  the  inlet  Nernst  potentials  should 
be  identical  if  the  inlet  fuel  composition  is  the  same  in  both 


X(cm) 

Fig.  3.  Plot  of  Nernst  potential  as  a  function  of  cell  axial  position  from  the  fuel  inlet  (i.e.  closed-end  of  cell)  at  low  average  current  density  (10 mA/cm2) 
and  high  average  current  density  (520  mA/cm2). 
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Fig.  4.  Experimentally  measured  fuel  utilization  sensitivity  curves  for  average  current  density  of  300mA/cm2  with  overlaid  numerical  fit  for  three  different 
leak  scenarios. 
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cases.  The  answer  lies  in  the  fact  that  the  same  number  of 
elements  (50)  was  used  for  both  the  high  average  current 
density  and  the  low  average  current  density  cases.  However, 
in  the  low  average  current  density  case,  the  fuel  is  utilized 
very  rapidly  along  the  length  of  the  cell.  The  first  data  points 
in  both  curves  do  not  represent  the  Nernst  potential  at  the 
absolute  closed-end  (which  is  the  fuel  inlet  point)  but  an  av¬ 
erage  across  the  first  element.  Since  the  rate  of  decrease  of 
Nernst  potential  near  the  inlet  is  much  more  rapid  for  the 
low  average  current  density  case,  the  average  Nernst  poten¬ 
tial  across  the  first  element  for  this  case  is  lower.  Indeed,  if 
the  calculation  is  repeated  with  a  larger  number  of  elements 
(200)  the  first  data  point  increases  from  920  to  960  mV  for 
the  low  average  current  density  case. 

Fig.  4  shows  the  fuel  utilization  sensitivity,  i.e.  the  first 
derivative  of  the  terminal  voltage  of  the  cell  with  respect  to 
fuel  utilization,  versus  the  fuel  utilization  for  three  different 
fuel  leak  rates  at  an  average  current  density  of  300  mA/cm2 
along  with  the  experimentally  measured  data.  As  can  be 
seen  from  Fig.  4,  the  overall  trend  of  these  curves  is  very 
similar  to  the  curves  presented  in  Fig.  2.  Over  the  course  of 
our  experiments  it  has  been  found  that  the  threshold  value 
of  fuel  utilization,  i.e.  the  fuel  utilization  at  which  a  signif¬ 
icant  change  in  the  slope  of  the  derivative  occurs,  is  very 
sensitive  to  the  fuel  leak  rate  of  a  given  cell;  higher  the  leak 
rate  the  lower  the  threshold  value  of  the  fuel  utilization.  In 
the  present  instance,  a  leak  rate  of  20  cm3 /min  fits  the  mea¬ 
sured  experimental  data  very  well.  Thus,  the  numerical  re¬ 
sults  presented  here  in  addition  to  corroborating  the  main 


details  of  the  simple  analytical  model  presented  earlier  can 
also  be  used  to  calculate  the  leak  rates  of  tubular  SOFCs 
through  measurement  of  the  sensitivity  of  the  terminal  volt¬ 
age  to  the  fuel  utilization. 

3.  Conclusions 

In  the  present  work,  we  have  provided  analytical  and  nu¬ 
merical  calculations  to  model  fuel  utilization  sensitivity  tests 
on  tubular  SOFCs.  The  analytical  model  has  been  used  to 
qualitatively  explain  the  shape  of  the  fuel  utilization  sensi¬ 
tivity  curves  and  the  numerical  model  to  quantify  fuel  leak 
in  tubular  SOFCs.  The  analysis  presented  here  provides  the 
capability  to  further  understand  how  leaks  affect  the  ultimate 
performance  of  tubular  SOFCs. 
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